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^ ; ABSTRACT 

^ ■ We determined C, N and a-element relative abundances in the gas surrounding six 

^f) ' QSOs at an average redshift of < 2: >~ 2.4, by studying six narrow associated absorp- 

Tlj- , tion systems in UVES high-resolution spectra. We found five systems with a metallicity 

' (measured by C/H) consistent or above the solar value. The ionization structure ob- 

I served in the associated systems is clearly different from that of the intervening ones, 

, indicating that the associated systems are influenced by the strong UV flux from the 

^ ' QSO. There is a possible correlation (anticorrelation) between [N/C] ([Si/C]) and 

I [C/H] of the studied associated systems, and [N/C] > when [C/H] > 0. We have 

r-j . compared these observational results with the predictions of a model simulating the 

C^h' joint evolution of QSOs and their spheroidal hosts. The agreement turns out to be 

JL , very good, in particular, the case envisaging massive haloes and high star-formation 

5^ ■ rates recovers both the correlation between [N/C] and [C/H] and the anticorrelation 

I for [Si/C] vs. [C/H]. Narrow associated absorption systems prove to be powerful trac- 

^ . ers of the chemical abundances in gas belonging to high redshift spheroidal galaxies. 

' The outflow of this same gas, triggered by the QSO feedback, is probably going to 

. ^ \ contribute to the early enrichment of the surrounding intergalactic medium. A larger 

. statistics, possibly increasing the number of ionisation stages, chemical elements and 

' the redshift range, would allow us to put firm constraints on detailed chemical evolu- 

I tion models of galaxies at high redshifts. 
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1 INTRODUCTION 

In this work, we want to address the star formation history 
and the evolution of massive early-type galaxies at high red- 
shifts by measuring in a reliable way the metallicity and the 
chemical abundances of gas belonging to host galaxies and 
environments of QSOs. 

Once considered rare and exotic objects, QSOs could 
instead represent a necessary phase in the evolution of mas- 
sive early-type galaxies. This interpretation is supported by 
several pieces of evidence: Massive Dark Objects (MDOs, 
generally interpreted as dormant black holes) with masses 

* Based on material collected with the European Southern Obser- 
vatory Very Large Telescope operated on Cerro Paranal (Chile). 
Proposal 116.A-0106, 65.O-0299, 67.A-0078 and 69.A-0204 
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in the range ~ 10^ - 3 x 10^ Mq are present in essentially all 
local gala xies with a sub stantial spheroidal component (see 
'Kormcrid v fc Gebhardtl lSoOl. for a review), on the other 
hand the host galaxies of low redshift powerful AGN (radio- 
loud and radio-quiet QSOs and radio galaxies) are, in all 
the studied cases, lu minous elliptical galaxies with L > L* 
JPunlop et alJl2003l) . 

The observational properties inferred for cluster and 
field elliptical galaxies up to redshift 2 ~ 1 imply a high 
uniformit y and synchronization in the galaxy fo rmation pro- 
cess fe.g. lEllis et anilOQTtlBernardi et alJll998l). The evolu- 
tion with redshift of their optical-IR colours llstanford et alJ 
[1998) is consistent with the passive evolution of an old stel- 
lar population formed a,t z > 2 — 3 and the measured positive 
[Mg/Fe] elemental ratio can be explained by a short and in- 
tense star formation b urst (e.g. IWorthev. Faber fc Gonzalea 
119921: lMatteuccilll994l) . 
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In the standard framework of the hierarchical evolu- 
tion of structures in a cold dark matter (CDM) universe 
large objects form by a sequence of mergers of smaller 
proto-galaxies. In particular, massive ellipticals are gener- 
ated at low redshifts {z < 2) from the merger of two 
large disk galaxies which formed sta rs at a constant mod- 
erate rate up to that moment (e.g. iBaugh. Cole fc FrenkI 
119961: iKauffmann fc Chariot ' 1998). In the merging event, 
the black holes (BHs) pre-existing in the progenitor galax- 
ies coalesce and a fraction of the cold gas is accreted 
by the new BH which activates as a QSO, the rest 
of the c old gas is transformed into stars in a sud- 
den burst (Wang S^Biemwan 1998; Kauffmann & Hach neld 
I2OOOI: fVolonteri. Haardt fc Madau.2003, : .Menci et al., 200^ 

A different prescription in the framework of the hi- 
erarchical scenario is the anti-hierarchical baryon collapse 
where the formation of stars and of the central BH takes 
place on sh orter time-scales within more massive dark mat- 
ter haloes llMon aco^ Salucci fc Danese||2000|: [G ranato et alJ 
I2OOII: I Archibald et alJ2002l:lGranato et al.l2004l) . Supernova 
heating and QSO feedback are the physical processes that 
reverse the order of formation of galaxies compared to that 
of DM haloes because they slow down star formation most 
effectively in shallow potential wells. In the more massive 
DM haloes star formation goes on rapidly causing at the 
same time the growth of the central BH which accretes the 
cold gas slowed down by the radiation drag. When the QSO 
activates, strong winds originate sweeping the interstellar 
medium and halting both the star formation and the BH 
growth. The time delay between the star formation onset 
and the peak of the QSO activity is again shorter for larger 
haloes. For the most massive galaxies (Afhaio ^ 10^^ -^q) 
virializing at 3 < z < 6, this time is < 1 Gyr, implying 
that the bulk of star formation may be completed before 
type la supernovae have the time to significantly enrich the 
interstellar medium with iron. A detailed analysis of the 
chemical evolution e xpected for this model is reported in 
iRomano et all tOO'j) . 

The two above described scenarii predict different chem- 
ical abundances, in particular at redshifts larger than ~ 2. 
The metallicity and the elemental abundances of high red- 
shift galaxies are hard to measure; on the other hand, high 
and intermediate resolution spectra of QSOs at redshifts as 
large as z ~ 6 can be easily obtained with the present in- 
strumentation. 

We studied associated narrow^ absorption lines exploit- 
ing high resolution, high signal-to- noise ratio spectra of 
2 < z < 3 QSOs obtained with the UVES spectrograph 
and a model for the photoionisation of the gas to derive 
chemical abundances in the QSO environments. Our results 
suggest that at these redshifts the gas associated with the 
QSO and with its host galaxy has already been enriched by 
the products of an intense star formation episode. 

Section 2 introduces the diagnostics that we used to 
determine the chemical abundances in QSO environments 
and reports previous results. In Section 3 we describe the 



^ The adjective "narrow" is used to distinguish this class of 
absorptions from the Broad Absorption Lines characterised by 
FWHM > 2000 km s"i and arising in gas ejected by the QSO at 
large velocities (see also Sections 2 and 7) 



selection criteria and the characteristics of our sample of as- 
sociated narrow absorption line systems; the adopted pho- 
toionisation model and the methodology are reported in Sec- 
tion 4. Section 5 is devoted to the description of our results, 
which are compared with model predictions in Section 6. In 
Section 7, we summarise the results on QSO chemical abun- 
dances obtained using other methods. We draw our conclu- 
sions in Section 8. 



2 CHEMICAL ABUNDANCES IN THE 

VICINITY OF QSOS MEASURED WITH 
ASSOCIATED NARROW ABSORPTION 
LINES 

The "narrow" absorption lines (NALs) are the most numer- 
ous in QSO spectra. A practical way to define them is that 
they have to be narrow enough that important UV doublets 
are not blended, i.e. to have FWHM < 200 to 300 km s"^ 

Here we are concerned only with the NALs falling 
within ±5000 km s~^ of the systemic redshift: the so called 
"associated" systems (AALs). In particular, since we are in- 
terested in determining the metal abundances and the physi- 
cal conditions in the gaseous environments close to QSOs, we 
would like to identify which of these systems are also "intrin- 
sic", that is physically associated with the QSO. The intrin- 
sic nature of individual AALs can be inferred from various 
indicators, for example (1) time variability in the absorp- 
tion lines requiring dense, compact regions and thus intense 
radiation fields near the QSO for photoionisation, (2) high 
space densities measured directly from excited-state fine- 
structure lines, (3) partial coverage of the emission source 
measured via resolved, optically thick lines with too-shallow 
absorption troughs, (4) spectropolarimetry that reveals an 
unabsorbed or less absorbed spectrum in polarized light, (5) 
smooth and relatively broad absorption line profiles that are 
unlike intervening absorbers, and (6) higher ionisation states 
than intervening absorbers. 

Previous studies on AALs have found that 
in general these syste ms have solar or supersolar 

metallicities (e.g. [^^^mglg^^g^gr on fc Pet ^jgrnl 

19931: iMaller. .Takobsen fc PerrvmanI Lli 



Pet^ear^^^u^l^fe^ CarswcU 1994; Hamann et al.' "1997 
Petitiean fc Srianand. 1999, : .Srianand fc Petitiean..200fl ') . In 



a few cases moderately subsolar metallicities are observed, 
but they are still significantly higher than those f ound in 
intervening systems at the same redshifts (Savagl Fo et alJ 
Il997fl . Besides, a marked change in the metallicity of QSO 
absorption systems from values smaller than 1/10 solar to 
solar or larger values is observed at a blues hift of ~ 15000 
km s"'^ relative to t he QSO omission lines ijPetitiean et alJ 
ll994l:lFranceschini fc Gratton 1992). 

In order to derive reliable metallicities and chemi- 
cal abundances we had to carefully select our AALs. The 
adopted criteria and the sample are described in the follow- 
ing section. 
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3 THE SAMPLE 

3.1 Selection criteria and data analysis 

We have considered the 16 QSO spectra of the ESO Large 
Programme "The Cosmic Evolution of the Intergalactic 
Medium" (116.A-0106A, P.I. J. Bergeron) available to the 
p ublic and the 6 QSO spectra of our pair archive described 
in lD'Odorico. Petitiean fc Cristianil ll2002h . All of them were 
obtained with the UVES spectrograph at the Kueyen unit 
of the ESO VLT (Cerro Paranal, Chile) with a resolution 
R ^ 40000 and a wavelength coverage 310 < A < 1000 nm. 
The signal to noise ratio (S/N) of the Large Programme 
spectra is S/N ~ 50, while in our spectra it is varying be- 
tween S/N ~ 5 and 20 depending on the brightness of the 
observed QSO. We used the Large Programme spectra re- 
duced with the automatised, refined version of the UVES 
pipe line (Balloster et al. 20001*) devised by Bastien Aracil 
(see lAracil et al .12004) . The normalisation has been carried 
out by manually selecting spectral regions not affected by 
evident absorption and by interpolating them with a spline 
function of 3rd degree. 

In these spectra we looked for C iv and N v doublets 
within ±5000 km s~^ from the emission redshift of the QSO. 
If both absorptions were detected at the same redshift we 
considered the system as a candidate intrinsic absorber and 
we looked for other associated ions (like C ll, C ill, N ill. Si ll, 
Sim, Si IV, Mgii, etc). 

Among the 22 QSOs in the sample, 16 showed asso- 
ciated C IV doublets summing up to a total of 34 systems 
of which 15 had detectable Nv lines too. In 11 of the 19 
systems without N v, the rest equivalent width of the C iv 
A1548 absorption hue was lower than ~ 0.03 A. Since in 
general the Nv/Civ equivalent width ratio is lower than 
one, the Nv transition even if present could be under our 
detection limit. Four of the 8 systems with stronger Civ 
doublets, have C lll/C iv equivalent width ratios larger than 
one. As we will discuss in Section 5 this is a characteristic of 
intervening systems while the contrary is true for associated 
ones. The remaining 4 systems do not show hints of intrin- 
sic nature by any of the six properties listed in the previous 
section, except one for which an effect of partial coverage is 
detected. 

In order to constrain the photoionisation model (see 
Section 4) and to get reliable results for the metallicity and 
the element relative abundances we required also that at 
least two ionisation states of the same chemical element were 
detected and that the corresponding absorption lines were 
not saturated. This requirement reduced the number of vi- 
able systems to 6. 

The systems excluded from the sample by our defin- 
ing criteria can be divided into two groups: the first one 
is formed by those systems for which we did not observe 
enough ionic transitions to constrain the photoionisation 
model, for the majority of them only a weak C iv doublet 
was detected. In some cases, more lines could be observed 
extending the spectral wavelength coverage with new obser- 
vations. The second group is made by a few strong systems 
where most of the lines are saturated not allowing the de- 
termination of reliable column densities. 

We fitted the observed absorption lines with Voigt pro- 
files in the LYMAN context of the MIDAS reduction pack- 
age (iFontana fc Balleste J [l995l) . Adopted atomic parame- 



Table 1. Atomic parameters used in the fitting of the fines which 
differ from Mort on 1 1991) 



Ion 


Rest wavel. 


Osc. strength 




(A) 




C IV. . . 


1548.204(2)1 


0.1908 




1550.781(2)1 


0.09522 


Mgii. . . 


2796.3543(2)2 


0.6123 




2803.5315(2)2 


0.3050 


Alii. . . 


1670.7886(1)1 


1.833 


Si 11. . . 


1304.3702 


0.086^ 




1526.70698(2)1 


0.1 10^ 




1808.01288(1)1 


0.0022-* 


Si IV. . . 


1393.76018(4)1 


0.5140 




1402.77291(4)1 


0.2553 



1 Griesmann & Kling (2000); ^ Pickering et al. (2000; 2002); 3 
Spitzer & Fitzpatrick (1993); * Bergeson & Lawler (1993) 

Table 2. Relevant parameters of the QSOs studied in this work. 
The absolute magnitudes are taken from the 10th edition of the 
QSO catalog by Veron-Cetty & Veron (2001) 



QSO Name 




Ma 


gnitude 


Abs. Mag. 


UM680 


2.1439 


V= 


■18.6 


-27.5 


UM681 


2.1219 


Bj 


= 18.8 


-27.3 


HE1158-1843 


2.453 


V= 


16.93 


-29.4 


Q2343+1232 


2.549 


V= 


17.0 


-29.7 


Q0453-423 


2.661 


v= 


d7.06 


-29.2 


PKS0329-255 


2.685 


v= 


d7.51 


-29.2 



ters differing from those reported in lMortonI lll99 J) are given 
in Table [H 

In the following we will describe in detail the selected 
systems. The properties of the background QSOs are re- 
ported in Tabled 

3.2 Individual systems 

System at z^bs = 2.123 towards UM680 

The system at Zabs ~ 2.123 (zcm — Zabs ^ 1990 km s~^) has 
strong Civ, Nv, Ovi and Siiv lines. Also Cii, Nil, Mgii, 
Alll, Sin and Siv are observed. The central components of 
C IV are too heavily saturated and for O vi the signal-to- 
noise ratio is too low to determine reliable column densities. 
Also the Hi Ly-a line is saturated. In order to obtain an 
estimate of the Hi column density for the velocity compo- 
nent corresponding to the metal absorptions, we fixed its 
redshift as the average of the two Si ll and Si iv component 
redshifts. The velocity profiles of the Ly-a and Ly-/3 tran- 
sitions are not constraining the fitting process which does 
not converge to a unique result. We had then to adopt a 
Doppler parameter for the central component, while leaving 
the other component parameters free. For a typical value 
6 = 20 km s~^ observed for Ly-a lines at redshift around 
2.1 (see iKim et al. 2001) a logiV(Hl) = 15.2 is derived. To 
estimate the uncertainty on this determination, we deter- 
mined the Hi column densities for two limit values of the 
Doppler parameter: 6 = 15 km s~i (consistent with the lower 
limit observed in the 6-distribution of the Lyman forest) and 
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-150 150-150 O 150 

Velocity (km/s) 

Figure 1. Ionic transition lines observed at redsliift z = 2.123 
(marked by the vertical dotted line) in the spectrum of QSO 
UM680. The result of the best fitting for the analysed absorptions 
is overplotted on the spectrum. The short-dashed line represents 
the noise 



Table 3. Parameters obtained fitting the absorption system at 
^abs ~ 2.123 in the spectrum of UM680 



Ion 


Redshift 


log N 
(cm~^) 


b 

(km s-i) 


Hi. . . 


2 


122327 ± 0.00005 


13.7 it 0.2 


30 it 3 




2 


12247 ± 0.00003 


1 Q 2 -1- n s 

-LvJ. _1_ W. vJ 


15 it 4 




2 


12302(1)1 


15.2^ 


zu^i; 




2 


123586 ± 0.000008 


1 Q no + n r)3 


19 ± 1 




2 


12399 ± 0.00004 


13.0 it 0.1 


21 ± 4 


C II. . . 


2 


12295(f) 


14 it 0.2 


4.5 it 0.6 




2 


1231(f) 


14.6 it 0.2 


6.2 it 0.6 


Nil. . . 


2 


123110 ± 0.000009 


13.79 it 0.06 


6.2 ±0.6 


Nv... 


2 


122445 ± 0.000008 


13.21 it 0.03 


16 ± 1 




2 


12293 ± 0.00002 


13.7 it 0.1 


4± 1 




2 


123058 ± 0.000008 


14.42 it 0.05 


9.1 ±0.5 




2 


123284 ± 0.000008 


12.77 it 0.07 


3(f) 




2 


12357 ± 0.00002 


12.93 it 0.07 


14 ±3 


Mgii. . . 


2 


122941 ± 0.000004 


12.51 it 0.02 


3.6 ±0.4 




2 


123091 ± 0.000004 


13.11 it 0.03 


5.2 ±0.3 


Alii. . . 


2 


12306 ± 0.00002 


11.25 it 0.09 


6.2 ±0.6 


Sin... 


2 


122949 ± 0.000005 


12.29 it 0.06 


2.5 ±0.9 




2 


123101 ± 0.000002 


13.03 it 0.02 


5.5 ±0.3 


Si IV... 


2 


122940 ± 0.000007 


13.18 it 0.06 


4 ±0.7 




2 


123091 ± 0.000007 


13.9 it 0.1 


6.4 ±0.6 


S IV. . . 


2 


123091 ± 0.000009 


14.39 it 0.05 


11 ± 1 



1 We indicate with an (f) those parameters that were assumed 
and not determined by the fitting process 

^ For the uncertainty on this determination see the discussion in 
the text 
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Figure 2. a) Theoretical Si II* (1264)/Si ll(1260) column density 
ratio as a function of temperature at increasing electron density 
(from bottom to top, as marked on the picture). The dashed lines 
represent the upper limits on the Si II* /Si II column density ratio 
(horizontal line) and on the temperature (vertical line) for the 
absorption system at Za^s = 2.123 in the spectrum of UM680; 
b) Theoretical C II* (1335)/C ll(1334) column density ratio as a 
function of temperature. Here the dashed lines correspond to the 
limits obtained for the absorption system at ^^bs — 

2.122 in the 

spectrum of UM681 



6 = 30 km (the maximum value allowed by the Ly-/3 ve- 
locity profile). They turned out to be logA'^(Hl) = 16 and 
log A'^(Hl) = 14.8 respectively. In the following computation 
of chemical abundance ratios for this velocity component we 
will adopt as the error on the Hi column density 0.6 dex to 
take into account this uncertainty on its determination. 

All the measured column densities are reported in Ta- 
ble]^ and the ionic transitions are shown in Fig. 
The fine-structure transition absorptions C ll* correspond- 
ing to the strong absorption lines due to C ll A 1036, 1334 
are absent. Besides, we did not observe the Sin* A 1264 fine- 
structure line associated with the Sill A 1260 absorption. 
The upper limits on the Sin* /Sill column density ratio for 
the two observed components are 0.19 and 0.04. The cor- 
responding upper limits on the temperature derived from 
the Doppler parame ters are log T < 4 and log T < 4.7 re- 
spectively. Following ISrianand fc PetitieanI i2000l) we deter- 
mined an upper limit on the electron density of rie ~ 10 — 30 
cm"'^ applying the average values of the previously re- 
ported limits to the formul a of pure coUisional excitation 
iFitzpatrick fc SpitzerlllQQTl . and references therein, note 
that excitation by hydrogen atoms is unimportant): 

N{X ir) ^ r^e 712(e) 

7V(XII) A21 + Tie 721(e)' ^ ' 

where X is either silicon or carbon, 7 are the collision ex- 
citation and de-excitation rates, and A21 is the radiative 
decay rate (see Fig. |5J . Since the gas is optically thin it is 
reasonable to assume ric « nu- 
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Figure 3. Ionic transition lines observed at redsliift z = 2.122 
(marked by the vertical dotted line) in the spectrum of QSO 
UM681. The result of the best fitting for the analysed absorptions 
is overplotted on the spectrum. The short-dashed line represents 
the noise 









II 1 'Ly y '. 


PTf -l-rTr-l'i- 






II I Ly a 






C IV iD-10 . 






1 1 i 1 1 1 1 

r IV I faTiO : 






cm 977 





-150 O 150 -150 O 150 

Vnloc'it.y (kT7-|/s) 

Figure 4. Ionic transitions lines observed at redshift z ~ 2.44264 
(marked by the vertical dotted line) in the spectrum of QSO 
HE1158-1843. The zero level of the Civ and Nv doublets has 
been modified to coincide with the minimum flux of the stronger 
line of the doublet in order to carry out the fitting. The same v^as 
done for the Ly-a absorption line. The result of the best fitting 
for the analysed absorptions is overplotted on the spectrum. The 
short-dashed line represents the noise 



An absorption system at Zabs = 2.122 (^om — Zabs ~ — 10 
km s~^) is detected with lines due to Civ, Nv, Ovi and 
Si IV. The low ionisation transitions of Cll, Nil, Sill, Alll 
and Mg ll are present too (see Fig. |3J . The H I Ly-a absorp- 
tion is saturated and the Ly-/3 line is partially blended. To 
measure the H l column density we fixed the redshift of the 
main component as the average of the two components ob- 
served in the low-ionisation transitions. We used only the 
Ly-a line to determine the fitting parameters and verified 
the consistency of the result with the blue wing of the Ly-/3 
velocity profile which looks free from major blending. All 
the measured column densities are reported in Table 2] 
Also for this absorption system we did not detect the fine- 
structure transition lines of C ll and Si ll. We obtained the 
upper limits on the column density ratio C ll* A 1335 over 
C II A 1334 and on the gas temperature for the two compo- 
nents of the system: C ll*/C ll < 0.16 and log T < 4.87, and 
C II* /C II < 0.18 and log T < 4^6 respectively. In Fig. HI we 
plot the curves of the C ll ratio as a function of temperature 
for constant ric as given by eq. Q The observed limits are 
consistent with an electron density of nc < 1 — 5 cm"''. 
Note that this absorption system and the one along the line 
of sight to UM680 are less than ~ 100 km apart in 
redshift and that the two lines of sight are separated by 
~ 700 kpc proper distance at this redshift. 

System at 2abs ~ 2.4426 towards HE1158-1843 

A ~ 350 km s~^ wide system is observed at an average ve- 
locity of ~ 900 km s~^ from the systemic redshift of QSO 
HE1158-1843 with associated C iii, C iv, N v and strong O vi 
absorption lines, as shown in Fig. 0] N ill. Si ill and S iii are 
blended if present, while Siiv and Siv are not observed. 
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Figure 5. Civ doublet observed at z ~ 2.4426 in the spectrum 
of QSO HE1158-1843. The thick solid line represents the profile 
fitting based on the intensity of the C IV A 1548 A absorption 
only. When applied to the companion line of the doublet it clearly 
shows the effect of partial coverage 
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Table 4. Parameters obtained fitting the absorption system at 
Zabs ~ 2.122 in the spectrum of UM681 



Ion 


Redshift 


log N 


b 






(cm~^) 


(km s-i) 


Hi... 




13 ± 0.3 


29 ± 7 




2 1 22(11 ^ff 1 


15.7 ± 0.2 


26 ± 1 




2 1 22500 + 000008 


14 16 + 02 


22 ± 1 




2.1229 ± 0.0001 


12.4 lb 0.4 


17 ± 9 


C 11. . . 




1 Q on _L fl 1 


11 ± 2 




9 1 991 1 n nnnno'^ 


1 Q 41 _|_ n 07 


6 ± 1 


Civ... 


2.12183 ±0.00002 


13.6 ±0.2 


21 ±3 




2.122016(f) 


15.23 ±0.04 


13.5(f) 




2.12235 ±0.00002 


13.62 ±0.06 


23 ±4 




2.122520 ± 0.000002 


13.3 ±0.1 


9±1 


Nv. . . 


2.1220193 ±0.000004 


13.49 ±0.03 


13 ± 1 


Ovi. . . 


2.122006 ± 0.000009 


14.5 ±0.03 


28 ± 2 




2.12254 ±0.00002 


13.92 ±0.06 


18 ± 0.2 


Mgii. . . 


2.121932 ± 0.000008 


12.1±0.1 


11 ±2 




2.122128 ± 0.000003 


12.38 ±0.07 


6±1 


Alii. . . 


2.121932(f) 


11.40 ±0.07 


11(f) 




2.122128(f) 


11.28 ±0.08 


6(f) 


Si 11. . . 


2.121930 ± 0.000008 


12.54 ±0.09 


11 ±2 




2.122113 ± 0.000002 


12.68 ±0.07 


6±1 


Si IV. . . 


2.122019 ± 0.000001 


14.12 ±0.05 


14.4 ±0.5 



In Fig. 13 we show that the intensity ratio of the Unes of 
the Civ doublet is not 2 as expected, this could be due to 
an incomplete coverage of the continuum source by the ab- 
sorber. The unnatural ratio of the N v doublet lines and the 
disagreement between the Hi Ly-a line and the Ly-/3, Ly-7 
line intensities confirm the hypothesis of partial coverage. In 
order to measure column densities for the observed ions we 
artificially fixed the zero level at the bottom of the stronger 
lines of the C iv and N v doublets and carried out the fitting. 
We followed the same procedure for the Ly-a for which we 
adopted the redshifts of the main C iv components. The col- 
umn density obtained for the central component of Ly-a is 
consistent with the velocity profiles of Ly-/3 and Ly-7. They 
are reported in Table |5] 

The upper limits on the column densities for the un- 
detected Si IV and Siv transitions are logA'^(Silv) < 11.35 
and log A'^(S iv) < 12.9 where we have adopted the Doppler 
parameter of C iv for the component at Zaba — 2.4426. 
The partial coverage effect suggests that the size of the ab- 
sorbing cloud is comparable with that of the continuum 
source, in the hypothesis of spherical geometry it is pos- 
sible to derive a lowe r limit to the absorber total density 
JPetitiean et alJll994ft : 

- > 5 " ('-^) " (2) 

\^ 1014-2 cm-2 y \^ 10-5 y \^PcJ 

where we have adopted rcont = 1 pc as the extension of the 
ionizing source considering that both the broad line region 
and the continuum emitting region contribute to the ionizing 
flux. 

We computed the chemical abundances only for the 
strongest component at Zabs ~ 2.4426 for which the column 
densities of both H I and metals where best determined. 

System at Zabs = 2.5697 towards Q2343+1232 

A system with three main components is observed in the 



Table 5. Parameters obtained fitting the absorption system at 
^abs ^ 2.4426 in the spectrum of HE1158-1843 



Ion Redshift log N b 









(cm 




(km s-i) 


2 


.44061 ± 0.00002 


12 


.02 


± 


0, 


,05 


20 (f) 


2 


.44195 ± 0.00002 


13 


.72 


± 


0, 


,04 


33 ± 2 


2 


.442263 (f) 


13, 


.65 


± 


0, 


,05 


23.2 ± 0.8 


2 


.442616 (f) 


14, 


.22 


± 


0, 


,02 


29.4 ± 0.3 


2 


.443143 (f) 


13 


.56 


± 


0, 


,01 


22.0 ± 0.3 


2 


.443506 (f) 


13, 


.38 


± 


0. 


,01 


21.8 ± 0.1 


2 


.444044 (f) 


12 


.87 


± 


0, 


,01 


41.0 ± 0.9 


2 


.441214 ± 0.000008 


13 


.22 


± 


0, 


,02 


31 ± 1 


2 


.44121 ± 0.00004 


11 


.70 


± 


0, 


,09 


15 ± 3 


2 


.441768 ± 0.000004 


12 


.24 


± 


0, 


,02 


6.9 ± 0.6 


2 


.44214 ± 0.00001 


12 


.28 


± 


0, 


,02 


13.2 ± 0.9 


2 


.442643 ± 0.000005 


12 


.43 


± 


0, 


,02 


6.8 ± 0.5 


2 


.44308 ± 0.00001 


12 


.38 


± 


0, 


,02 


13.7 ± 0.9 


2 


.443539 ± 0.000004 


12 


.26 


± 


0, 


,02 


8.2 ± 0.4 


2 


.44075 ± 0.00001 


12 


.55 


± 


0, 


,03 


20 ± 1 


2 


.441137 ± 0.000006 


13, 


.12 


± 


0, 


,02 


12.6 ± 0.4 


2 


.441411 ± 0.000006 


12 


.90 


± 


0, 


,04 


15 ± 1 


2 


Ainu ± 0.000001 


13 


.31 


± 


0, 


,01 


10.6 ± 0.3 


2 


.44211 ± 0.000006 


13 


.80 


± 


0, 


,03 


17.2 ± 0.6 


2 


AiHb ± 0.000006 


13 


.44 


± 


0, 


,06 


11.2 ± 0.4 


2 


.442627 ± 0.000006 


14, 


.05 


± 


0, 


,05 


10.7 ± 0.1 


2 


A ,10000 1 00000/^ 

.442898 ± 0.000006 


13 


.48 


± 


0, 


,03 


10.5 ± 0.6 


2 


.443131 ± 0.000006 


13 


.57 


± 


0, 


,01 


9.7 ± 0.2 


2 


A A A r\f\ 1 OOOOO"! 

.443499 ± 0.000001 


13 


.56 


± 


0, 


,03 


12.3 ± 0.1 


2 


.444045 ± 0.000002 


12 


.80 


± 


0, 


,01 


17.7 ± 0.6 


2 


.444378 (f) 


12 


.15 


± 


0, 


,02 


10 (f) 


z 






.05 


± 


0, 


,04 


5 (t) 


2 


.440667 ± 0.000008 


12, 


.58 


± 


0, 


,03 


11 ±1 


2 


.441144 ±0.000002 


12, 


.98 


± 


0, 


,04 


6.8 ±0.5 


2 


.441194 ±0.000008 


13, 


.54 


± 


0, 


,01 


25.5 ±0.9 


2 


.441780 ± 0.000008 


13, 


.44 


± 


0, 


,03 


13.0 ±0.5 


2 


.442101 ± 0.000008 


14, 


.03 


± 


0, 


,02 


16.8 ±0.6 


2 


.442291 ± 0.000002 


13, 


.73 


± 


0, 


,03 


8.9 ±0.3 


2 


.4426484 ± 0.0000008 


14, 


.33 


± 


0, 


,05 


11.1±0.1 


2 


.442886 ± 0.000008 


14, 


.31 


± 


0, 


,01 


30.5 ±0.5 


2 


.443153 ± 0.000001 


13, 


.42 


± 


0, 


,01 


3.5 ±0.2 


2 


.443519 ± 0.000001 


13, 


.82 


± 


0, 


,01 


10.0 ±0.1 


2 


.443752 ± 0.000008 


12, 


.51 


± 


0, 


,04 


4.2 ±0.6 


2 


.444066 ± 0.000003 


13, 


.08 


± 


0, 


,01 


14.7 ±0.7 


2 


.444396 ± 0.000008 


12, 


.57 


± 


0, 


,03 


10 (f) 


2 


.444680 (f) 


12, 


.20 


± 


0, 


,03 


8 (f) 



spectrum of Q2343+1232, redshifted by ~ 1700 km s"^ with 
respect to the QSO emission lines. It shows absorptions due 
to C III, C IV, N V and a weak Si iv. N lll is detected only at 
the redshift of the strongest component. The upper limits 
on the other two main components obtained adopting the 
Doppler parameter of C ill are: log A'^(N ill) < 13 at Zabs ~ 
2.568876 and logAr(Niii) < 13.1 at Zabs = 2.569432. The 
H I column density was obtained by the simultaneous fitting 
of the Ly-a and Ly-/3 absorption complexes which look free 
from blending. The detected ionic transitions are shown in 
Fig. |S] and the measured column densities are reported in 
Table il 

System at Zabs = 2.63618 towards Q0453-423 

A Civ doublet is detected at this redshift (zcm — 2abs~ 2040 
km s"^) showing a very weak associated Nv A 1238 ab- 
sorption and also weak Sim and Siiv (see Fig. |7|l. Ovi 
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Figure 6. Ionic transition lines observed at redshift z = 2.5695 
(marked by the vertical dotted line) in the spectrum of QSO 
Q2343+1232. The result of the best fitting for the analysed ab- 
sorptions is overplotted on the spectrum. The short-dashed line 
represents the noise 



Table 6. Parameters obtained fitting the absorption system at 
^abs =2.5697 in the spectrum of Q2343+1232 



Ion 


Redshift 


log N 
(cm~-^) 


b 

(km s^-"-) 


Hi. . . 


2 


568871(f) 


13.42 ± 


01 


17.7 ±0.5 




2 


569381 ± 0.000008 


13.65 ±0 


02 


22 ± 1 




2 


569957 ± 0.000008 


13.99 ±0 


01 


19.5 ± 0.6 




2 


570681 ± 0.000008 


13.34 ±0 


02 


24 ± 1 


Cm. . . 


2 


568868 ± 0.000008 


12.89 ±0 


04 


10 ± 1 




2 


56947 ± 0.000008 


13.15 ±0 


03 


14 ± 1 




2 


569936 ± 0.000004 


13.70 ±0 


05 


9(f) 


C IV. . . 


2 


568844 ± 0.000003 


13.19 ±0 


05 


4.2 ±0.5 




2 


568909 ± 0.000006 


13.31 ±0 


04 


11.6 ± 0.5 




2 


569232 ± 0.000006 


12.36 ±0 


06 


3(f) 




2 


569424 ± 0.000005 


13.43 ± 


03 


7.6 ±0.6 




2 


569587 ± 0.000006 


12.99 ±0 


07 


5.4 ±0.8 




2 


56990 ± 0.00001 


13.78 ±0 


08 


9.9 ±0.6 




2 


570007 ± 0.000006 


13.90 ±0 


07 


6.0 ±0.4 


Niii. . . 


2 


56994 ± 0.00002 


13.46 ± 


05 


13 ± 2 


Nv. . . 


2 


568869(f) 


13.20 ±0 


05 


6.9 ±0.8 




2 


568934(f) 


13.28 ±0 


05 


15 ± 1 




2 


569434 ± 0.000008 


13.34 ±0 


02 


14 ± 1 




2 


56993 ± 0.00002 


13.75 ±0 


09 


9.2 ±0.6 




2 


570023 ± 0.000008 


13.46 ± 


16 


4±2 


Si IV... 


2 


568852 ± 0.000003 


11.73 ±0 


09 


4.2 ±0.5 




2 


569432 ± 0.000005 


11.8 ±0 


1 


7.6 ±0.6 




2 


569596 ± 0.000007 


11.6 ±0 


1 


5.4 ±0.8 




2 


56990 ± 0.00001 


12.29 ±0 


08 


11 ± 2 




2 


570019 ± 0.000007 


12.24 ±0 


07 


4.8 ±0.8 
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Figure 7. Ionic transition lines of the absorption system at red- 
shift z = 2.63618 (marked by the vertical dotted line) in the 
spectrum of QSO Q0453-423. The result of the best fitting for the 
analysed absorptions is overplotted on the spectrum. The short 
dashed line represents the noise 



Table 7. Parameters obtained fitting the absorption system at 
2abs =2.63618 in the spectrum of Q0453-423 



Ion 


Redshift 


log N 


b 






(cm~^) 


(km s^-"-) 


Hi. . . 


2.6361797 ± 0.0000008 


14.76 ± 0.01 


20.0 ±0.1 


Civ... 


2.636071 ± 0.00002 


12.7 ±0.1 


12.4 ±0.8 




2.636187 ± 0.000001 


13.60 ±0.01 


8.7 ±0.1 




2.636484 ± 0.000006 


11.60 ±0.04 


3(f) 


Nv... 


2.636167 ± 0.000008 


12.43 ± 0.03 


12 ± 1 


Sim... 


2.636205 ± 0.000003 


11.72 ± 0.02 


7.6 ±0.4 


Si IV... 


2.636187 ± 0.000001 


12.34 ±0.01 


7.9 ±0.2 



and C III are possibly present although they are blended. 
The upper limit on the N ill and C III column densities 
adopting the redshift and Doppler parameter of Sim are 
logAr(Nlll) < 12.8 and logAf(Clll) < 13.55. The Hi Ly-a 
and Ly-/3 lines were used to determine the column density: 
logAf(Hl) ~ 14.763 ± 0.004 which is consistent with the 
blended Ly-7 absorption. All the measured column densi- 
ties are reported in Table Q 

System at Zabs = 2.708968 towards PKS0329-255 

This system is redshifted by ~ 1940 km s^^ from the QSO 
emission lines. It shows optically thin absorptions due to 
C IV, C III and N v. O vi is also detected, it is saturated 
and with at least two other components besides the one 
in common with the other ions (see Fig. |SJ|. The measured 
column densities for the observed transitions are reported 
in Table |H| The upper limits on the column densities of 
Niii and Siiv obtained adopting the Doppler parameter 
of Cm and Civ, respectively are: logAf(Nlll) < 12.5 and 
log7V(Silv) < 11.5. The Hi column density corresponding 
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Figure 8. Ionic transition lines of the absorption system at red- 
shift z = 2.70897 (marked by the vertical dotted line) in the spec- 
trum of QSO PKS0329-255. The result of the best fitting for the 
analysed absorptions is overplotted on the spectrum. The short 
dashed line represents the noise 



Table 8. Parameters obtained fitting the absorption system at 
^abs =2.7089 in the spectrum of PKS0329-255 



Ion 


Redshift 


log N 


b 






(cm~^) 


(km s-i) 


Hi. . . 


2.708917 ± 0.000008 


14.01 ±0.05 


15.9 ±0.8 




2.709304 ± 0.000008 


14.75 ±0.02 


32.6 ±0.6 




2.7096 ±0.0001 


13.6 ±0.15 


53 ±3.0 


Cm. . . 


2.70891 ± 0.00001 


12.07 ±0.05 


4±2 


C IV. . . 


2.708939± 


13.589 ± 0.005 


6.0 ±0.1 




2.709371 ± 0.000006 


12.47 ±0.03 


15 ± 1 


Nv. . . 


2.708913 ± 0.000002 


13.29 ±0.01 


5.5 ±0.2 


Ovi. . . 


2.70751 ± 0.00001 


12.81 ±0.04 


13 ±2 




2.70793 ± 0.00001 


13.09 ±0.08 


10 ± 1 




2.70824 ±0.00007 


12.9 ±0.1 


23 ±8 




2.70867 ±0.00001 


13.25 ±0.08 


7± 1 




2.708977 ±0.000003 


15.1 ±0.1 


9.5 ±0.8 




2.709390 ± 0.000004 


14.66 ±0.01 


15.6 ±0.6 




2.709800 ± 0.000004 


14.12 ±0.04 


8.0 ±0.5 




2.70998 ± 0.00002 


13.79 ±0.07 


14 ± 1 



to the metal absorption is well constrained by the simulta- 
neous fitting of the Ly-a, Ly-/3 and Ly-7 transitions which 
gives: logiV(Hl) ~ 14.01 ± 0.05. 



4 THE PHOTOIONISATION MODEL 

In order to obtain elemental abundances from the observed 
ionic column densities it is necessary to compute ionisa- 
tion correction s. To this purpose we used the code Cloudy 
l|Ferlandll200^ adopting as an ionising spec t rum a typical 
QSO spectrum derived from lCristiani fc Viol l)l99Cll and ex- 
trapolated in the region shortward of the Ly-a emission 



with a power law /(u) oc v'"'^ corresponding to the con- 
tinuum slope observed redward of the Ly-a emission. At 
energies higher than the Lyman l imit, we a. d opted a power 
law /(i/) (X following H amann et al.1 ^00^. This is 

consistent with the best observations of luminous QSOs and 
with the slope of the UV background at high energies. We 
verified that results did not change using a slightly differ- 
ent composite spectrum obtained from the FIRST bright 
QSO survey bv lBrotherton et al.l fcOOlT) . The absorption in 
the range 912 < Arest < 1216 A was corrected adopting 
the power law fitting the spectrum at larger wavelengths: 
/(i/) cx and for A.ost < 912 A we used oc 

as before. 

We assumed that components of different transitions 
which are at the same velocity arise in the same gas and we 
tried to reproduce with Cloudy all the ionic column densities 
observed at the same redshift with a single-region model. 

As a zero-order approach we computed for each system 
a grid of Cloudy models for a cloud with the measured Hi 
column density, solar metallicity and relative abundances, 
varying the ionisation parameter, which is defined as: 



where Li, is the luminosity density of the QSO spectrum, c 
is the speed of light, ro is the separation between the QSO 
and the illuminated face of the cloud and nu (cm~^) is the 
total hydrogen density (Hll + Hi + H2). 

In every studied system it was not possible to reproduce 
all the observed ionic column densities with a single U value, 
even allowing for a vertical shift to account for a non-solar 
metallicity. This was true also when considering groups of 
ions of similar ionisation state. Then, this discrepancy can- 
not be ascribed (only) to the fact that we are treating as a 
single region a multi-phase gas, and to be solved it requires 
the introduction of non-solar relative abundances. 

To derive the relative abundances and reproduce the ob- 
served set of column densities we followed three main steps. 

(1) The value of the ionisation parameter for the considered 
absorption system. Us, was fixed on the basis of the column 
density ratio of two, possibly close, ionisation states of one 
or more elements. 

(2) The column densities of all observed ions were then com- 
puted with Cloudy for U — Us, with the measured Hi col- 
umn density, solar metallicity and relative abundances. The 
ratios of the observed to the computed column densities were 
taken as the variations of the element abundances with re- 
spect to the solar values. 

(3) As a final step, we introduced new elemental abundances 
in the model computed from the variations determined in 
step (2). Then, instead of fixing Us we normalised the spec- 
trum with the absolute magnitude of the studied QSO and 
ran a grid of models varying the value of nn and finding the 
corresponding values of ro which gave viable solutions. 

Due to the degeneracy between the total density and the 
radius (refer to eq. 3), for all the reasonable values of wh it 
is possible to find a corresponding ro at which the set of col- 
umn densities matches the observed one. For three systems 
in Section 3.2 we derived physical limits on nn in an inde- 
pendent way which consequently gave limits on the distance 
from the source. In principle it is possible to derive a lower 
limit on nu, that is an upper limit on the distance from the 
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QSO, computing the total density of a cloud with the ob- 
served neutral hydrogen column density and a characteristi c 
dimension equal to the local Jeans length (see ISchavel200ll) . 
However, this approach, due to the low Hi column densities 
observed, does not provide stringent limits (ro,max ~ 10 — 20 
Mpc). 

We could not determine the metallicity of the gas from 
the abundance of iron, since we did not observe iron in our 
systems neither as Fell nor as Felll. Observations of stars 
in our Galaxy show th at [C/Fe1 is consistent with solar a t 
least for [Fe/H] > -1 JCarretta. Gratton fc Snedenll200Cl) . 
This is predicted by chemical evolution models of our Galaxy 
where a consistent amou nt of carbon is produced by low and 
intermediate mass stars jTimmes^Woo slev fc Weaveil[l995t 
IChiappini. Romano, fc Matteu^cil200^ enriching the inter- 
stellar medium on a time-scale similar to that of iron enrich- 
ment from Type la SNe. In our analysis, we used carbon as 
a proxy of iron and we computed the a-element abundances 
compared to it. 



5 RESULTS 

We investigated the star formation history in the close neigh- 
bourhood of six QSOs at emission redshifts between 2.1 and 
2.7 deriving the metallicities and relative chemical abun- 
dances of six associated narrow absorption systems. The re- 
sults of our calculations are summarised in Table 9 where 
errors on the abundance ratios are due mainly to the uncer- 
tainties in the column density determinations. 

1. For the systems associated with the QSOs UM680 and 
UM681 we were able to put upper limits on the total hydro- 
gen density in an independent way (see Section 3.2). They 
correspond in the Cloudy photoionisation models to lower 
limits on the distance from the emitting source of ro ~ 70- 
120 kpc and ro ~ 120-260 kpc for UM680 and UM681 re- 
spectively. The two absorption systems fall at very close red- 
shifts and present similar velocity structures, evidences that 
suggested the presence of a dif fuse gaseous structure pos- 
sibly including the two QSOs ISh aver fc RobertsonI 119831: 
[P'Odorico et al. 2002) . The metallicity of the two systems 
is quite uncertain due to saturated H I Ly-a lines. However, 
the system associated with UM680 has a [C/H] abundance 
ratio at least 5 times solar and larger than that of UM681 
which is consistent with solar. The [a/C] is larger in UM681 
than in UM680, while the two systems have [N/C] abun- 
dances consistent with solar. We speculate that the observed 
gas is an outflow of UM680 which is pierced by the UM681 
line of sight in an external, less enriched region. Deep imag- 
ing of the fleld could possibly shed some light on the nature 
of these absorbers. For example by detecting Ly-a-emitters 
at the same redshift of the absorbers tracing a large scale 
structure between the two lines of sight. 

2. For the absorber towards QSO HE1158-1843 it was possi- 
ble to compute a lower limit on the total density which trans- 
lates into a separation from the continuum source smaller 
than ~ 40 kpc. In this case we can assert that we are prob- 
ing the abundances in the interstellar medium of the host 
galaxy. The metallicity is about twice solar and the N/C 
ratio is slightly supersolar. Unfortunately we did not detect 
any a-element. 



Q0453-433 PKS03a9-S55 




Velocil.y (km/s) 

Figure 9. Observed intervening systems at Xabs = 2.4426 in the 
spectrum of Q0453-423 (left panels) and at z^bs = 2.4552 in the 
spectrum of PKS0329-255 (right panels) 

3. The three components of the system in QSO Q2343-f 1232 
show very similar relative abundances although the H I col- 
umn densities vary of a factor of 3 among them. The metal- 
licity and the ratio N/C are between 2 and 3 times solar and 
there is an indication of a-enhancement. 

4. The two systems at larger redshifts are characterised by 
small C IV equivalent widths and show different abundance 
patterns compared to the other systems. In particular they 
have undersolar N/C ratios. 

5. No significant correlation is observed between the veloc- 
ity separation from the quasar and the metallicity of the 
system. This suggests that the observed separations cannot 
be trivially related to the actual spatial distances between 
the central source and the absorber due both to uncertain- 
ties in the emission redshift determination and to peculiar 
velocities of the absorbing material. 

As shown in Fig. 1101 only one among the six stud- 
ied systems has a metallicity significantly lower than solar, 
Z ^ 1/6 Zq. The other 5 systems show values compara- 
ble or larger than solar. We conflrmed the supersolar N/C 
abundance ratio in those s ystems with Z > Zf?^ , as a lready 
found in other AALs (see iHamann fc Ferlandlll999l . for a 
review). On the other hand, in our AALs we measured en- 
hanced a-element/C abundance ratios at variance with the 
tentative detection of supersolar Fe/Mg abundance ratios in 
broad emission line regions (see Section 7) but in agreement 
with abundances measured in elliptical galaxies. This sug- 
gests that we are sampling regions where SNe la did not yet 
have the time to enrich the gaseous medium, implying that 
the bulk of star formation started less than ^ 1 Gyr before. 

Due to the lack of constraints on the total density of the 
absorbers we cannot reliably establish how close they are to 
the continuum source in all but one case. However, adopting 
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a density typical of the interstellar medium, wh = 10 cm~ , 
we obtained separations of the order of 100 kpc. To reach 
such a distance in 0.5 — 1 x 10* yrs, the average lifetime 
of a QSO, the gas should have travelled at a velocity of 
~ 1000 — 2000 km s~^, which is indeed of the order of the 
Av measured for our systems. 

Another way to test the hypothesis that these systems 
are under the strong influence of the associated QSO and 
probe chemical abundances in the QSO vicinity is to com- 
pare them with intervening systems with similar C iv equiv- 
alent width. In Fig. |^ we show two examples of intevening 
systems at velocity separations of ^ 18400 and 19300 km s^^ 
from the QSOs Q0453-423 and PKS0329-255 respectively. It 
is apparent that the ionisation pattern is completely differ- 
ent from the one observed in AALs, in particular the C ill 
transition is much stronger than the C iv A 1548 line while 
in the studied AALs the Clll/Civ column density ratio is 
always lower than one. This suggests that most optically 
thin intervening systems are ionised by the diffuse UV back- 
ground while associated systems receive the direct emission 
from the QSO. In particular, the ionisation energy of the 
C III — > C IV transition is close to a break present in the UV 
background spectrum at the He ionisation energy which is 
absent in the unabsorbed QSO spectra. 

The analysis of abundances in intervening systems is 
a difficult task because in many cases it is not possible 
to derive with confidence the Hi column density and, in 
general, there is an ambiguity on which kind of sources is 
ionising the gas. However, the vast majority of the studies 
on intervening Lyman limit and optically thin metal ab- 
sorption systems carried out up to now found undersolar 
metalhcities (e.g. iBergero n. et all[l994l: iKo hler et alJll99 9l: 
iProchaska fc Burleslll999l : lD'Odorico fc Petitieanll200ll) . Fi- 
nally, the damped Ly-a systems for which the determination 
of abundances is very precise due to negligible ionisation cor- 
rections, are characterised by an average metallicity lower 
than 1/10 solar at 2; > 1. 

We conclude that the associated systems show indeed 
a peculiar chemical abundance pattern and in the following 
section we try to frame our results into a model for the 
chemical evolution of large elliptical galaxies. 
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Figure 10. [N/C] and [Si/C] abundance ratios vs. [C/H] ob- 
tained for the analysed associated absorption systems and for the 
considered chemical evolution model. Solid triangles with error 
bars are our data. The abundances of the three components of 
the system towards Q2343-I-1232 have been averaged to a single 
value. The empty triangle in the bottom panel is the [O/C] abun- 
dance ratio for the system towards PKS0329-255 . Upper panel: 
solid d ots are the predictions of case B based on ICranato et aU 
i20n4ft for a DM halo mass of 3 X 10^^ Mq, zqsq = 2.1 and 
Zvir = 2.4, 2.5, 2.7 and 3 in order of increasing [C/H] abundance. 
The crosses are the predictions for a DM mass of 5 X 10" Mq, 
ZQSo = 2.55 and ^vir = 2.8 , 2.9, 3, 3.2, 3.5 . Solid squares are the 
results of case A based on iRomano et alJ j20o3) for a DM halo 
of Mjjaio = 1-37 X 10^^ '^O- Increasing sizes represent virializa- 
tion redshifts ^vir = 2.5, 3 and 4, the different groups of squares 
are the abundances at time-steps of 0.1 Gyrs from z^ir to zqso- 
Lower panel: Solid dots, crosses and solid squares are the same as 
above . Empty squares represent the predictions bv lRomano et all 
J2002l) for the abundance of [O/C] starting at z^ir = 3 



6 COMPARISON WITH MODEL 
PREDICTIONS 

We compared the present results on abundances with the- 
oretical predictions based on a physical mo del for the co- 
evolution of QSO and ho st galaxy systems llRomano et alJ 
l2002l: iGranato et al.l^2004^ . In order to illustrate the effects 
of the star-formation rate (SFR) history on abundances, we 
considered two cases for the gas distribution inside the viri- 
alized DM halo. The first one, case A, assumes that the gas 
closely follows the DM profile with no dumpiness, while in 
the second case, B, we assumed that after virilization on the 
average the gas follows the DM profile, but we introduced a 
clumping factor. As expected, in case A it is difficult to get 
a rapid star formation in DM haloes with Afh aio > 1.5 x 10^^ 
Mf7^ before the QSO shines (see Fig. 10 in iRomano et 
I2OO3), while in case B we can form stars very rapidly even 
in larger haloes (see lGranato et al]l2004^ . 

Both cases consider a single-zone galaxy and a double 



power-law IMF, i.e. $(M) oc M"°" for M < 1 Mq and 
$(M) (X M"^'^^ for 1 M0 < M < 100 Mq (in a notation 
where the Salpeter index would be 1.35). We also used for 
both c ases the chemical yields adopted bv IChiappini et al] 
i2003l their model 7). This yield set was chosen because it 
proved to give the best agreement with CNO osbervations 
for the Milky Way, the MlOl spiral galaxy, DLAs and dwarf 
irregular galaxies. 

For case A we simulate d a halo of mass Mhaio = 
1.37 X 10^^ Mq (model Id in IRomano et 3111200211 . We fol- 
lowed the evolution of gas abundances from the beginning of 
star formation to the shining of the QSO at time-steps of 0.1 
Gyrs and for different virialization redshifts. The model pre- 
dictions reported in Fig. llOl reproduce the correlation [N/C] 
vs. [C/H] and the a-enhancement, but the relation between 
[Si/C] and [C/H] is almost flat. 

The predictions for case B (see ICranato et al.ll2004h 
refer to the chemical evolution of two DM haloes of masses 
Mhaio = 3 X 10^^ and 5 x 10^^ Mq corresponding roughly 
to the absolute magnitudes Mb = —27.5 and —29.5 in the 
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hypothesis that the QSOs are shininR at the Ed dinRton lu- 
minosity and applying eq. 6 bv lFerraresd i2002l) : 



Mr 



108 Mp 



0.10 







1012 Mq 



(4) 



We explored virialization redshifts at time intervals of 
~ 0.5 and 1 Gyr from the QSO emission redhifts in our 
sample. As shown in Fig. 1101 the predictions cover the range 
of [C/H] where most of the observations lie. They reproduce 
both the observed correlation between [N/C] and [C/H] and 
the anticorrelation in [Si/C] vs. [C/H]. 

The differences between the predictions of case A and B 
are small. They are mostly due to the higher SFR that can 
be attained in case B before the QSO shining, due to the 
clumping factor that shortens the cooling time of the gas. It 
should be noticed also that this is a single-zone model, i.e. 
the results are averaged over the whole physical dimension 
of the galaxy, while a metallicity gradient is observed in 
elliptical galaxies. 

In conclusion, the high level of chemical enrichment and 
the a-enhancement observed in the QSO environments indi- 
cates that the massive elliptical galaxies hosting QSOs must 
have formed the bulk of their stellar population on short 
time-scales at high redshifts. 



7 ALTERNATIVE MEASUREMENTS OF 
CHEMICAL ABUNDANCES IN THE 
VICINITY OF QSOS 

Our results are in general agreement with QSO chemical 
abundances determined with other methods. 

Broad emission lines (BELs) observed in QSO spectra 
are t he m ost commonly used diagnostics to this purpose (see 
e.g. lOavidson fc Nctzcr 1972,; Hamann fc Fcrland 1999). 
BELs are known to originate in photoionised gas within 
~ 1 pc of the central continuum source. It has become the 
norm in BEL studies to take the N/O abundance ratio as 
a tracer of 0/H, that is of Z. The prominent metal lines, 
such as Civ A 1549, relative to Ly-a are not sensitive to 
the o verall metallicity for Z > 0.1 Zq iHamann fc Ferlandl 
[^Q^). On the other hand, observations in Hll regions indi- 
cate that the N/O oc O/H relation is valid for metallicities 
above ^ 1/3 to ~ 1/2 solar l| Shields"l976; Pagcl fc Edmundj 
198ll:lvan Zee. Salzer fc Havnes 19 98: Izotov fc Thuan 1999i: 
Pettiri^t^Lll2002[) . This abundance behaviour is attributed 
to "secondary" N production, whereby N is synthesized from 
existing C and O via CNO burning in stars. 

Calculations of BEL metallicities in large samples 
of QSOs spanning the redshift range < Zem < 5 
found typically solar or supersolar metallicities across 
the entire redshift range and no evidence of a de- 
crease at the highes t redshifts ( Ha mann fc Ferlarid| |l993^ 
Hamann et alJ [2OO2I: IWarner et al] l2002| iDietrich et alJ 
2003al : I Warner. Hamann fc DietrichI 12003). Precise esti- 
mates are difficult because metallicities derived from the N v 
lines (most notably N v/Hell) are typically ~ 30 % to a fac- 
tor of ~ 2 larger than estimates from the intercombination 
ratios (e.g. Niii]/Cni]). The reason for this discrepancy is 
not clear and also the absolute uncertainties are not easily 
quantified because they depend on the assumed theoretical 
models. 



The determination of the relative abundance of iron 
versus a-elements, which is taken as an indicator of the 
time elapsed from the beginning of the last star formation 
episode, relies for QSO BELs on the ratio Fen(UV)/Mgii 
A 2798, where "Fell(UV)" indicates a broad blend of many 
Fell lines between roughly 2000 and 3000 A which is very 
hard to measure. The tentative measurements carried out 
at low and interrnediato redshifts indicate Fe/Mg a factor of 
~ 3 above solar ([Wills et al. 1985) suggesting that SN la al- 
ready contributed to the gas enrichment. Recently, new mea- 
surements of the Fell/Mgll ratio were obtained for high and 
very high (z ~ 6) redshift QSOs showing no clear evolution 
with ti me (iFreudling. Corbin fc Koristal20o3 : IDietrich et alJ 
l2003hl) . However the consequences of these results on the 
early star formation history will be clearer only when the 
theoretical relationship between the observed Feii/Mgii 
emission ratio and the Fe/Mg abundance will b e assessed 
JVerner et alJll99gl l2003l : ISigut fc Pradhanll2003D . 

Also "broad" absorption lines (BALs) have been used 
to investigate the chemical and physical properties of the gas 
associated with QSOs. Indeed, BALs are believed to arise in 
material ejected by the QSO but still located very close to 
the central regions. The drawback in the use of these sys- 
tems is that broad profiles blend together all the important 
doublets and do not allow a reliable estimate of the column 
densities, in particular in the case of partial coverage of the 
source. However, the numerous studies on BALs indicate 



metallicities near or above the solar value (e.g.'lCorist^^^^ 
n996; H amann 1 99 ^:lArav et al.. 2001: Srianand fc Petitiean 
120011: iCupta et alJl2003l) ' 



8 CONCLUSIONS 

Up to now the main approach to study the chemical abun- 
dances in QSO environments has been the analysis of BELs 
observed in their spectra. Metallicities determined from 
BELs are consistent with solar or slightly supersolar values 
without a significant evolution in redshift. Other elemental 
abundances are very difficult to measure, in particular de- 
terminations of the ratio a/Fe are very uncertain. 

Associated narrow absorptions are complementary 
probes of the physical status of QSO-elliptical systems with 
respect to BELs. In general, they can be due to gas be- 
longing to the interstellar medium of the galaxy, outflowing 
under the effect of the QSO or re-infalling on the QSO itself. 
Furthermore, it is more straightforward to derive chemical 
abundances from absorption lines than from emission lines. 
We need only to determine and apply the proper ionisation 
corrections to convert the measured ionic column densities 
into relative abundances. 

In this paper, we selected six narrow absorption sys- 
tems lying within 5000 km s~^ from a QSO emission redshift 
and determined the abundances of C, N and a-elements in 
the gas they originate from. We used high resolution, high 
signal-to-noise UVES QSO spectra and applied a procedure 
based on the photoionisation code Cloudy to compute the 
chemical abundances starting from the measured column 
densities. 

We found that all systems but one in our sample have 
metallicities (measured by carbon) consistent with or larger 
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than solar. We found also a possible correlation of [N/C] 
and an anticorrelation of [Si/C] with [C/H] with supersolar 
values of [Si/C]. These results are suggestive of rapid en- 
richment due to a short star formation burst, of duration 
iburst ~ 1 Gyr (see Section 5). Since the very high lumi- 
nosity QSOs in our sample sh ould have Mbh > 10^ M g^ , 
assuming M^ph/MBH ~ 1000 iMcLure fc Dunlodl20o3) we 
expect SFR> 1000 Mq yr"^ in their hosts. 

The predictions of the model of chemical evolution for 
a spheroidal galaxy where the star formation depends on 
stellar and QSO feedback are in good agreement with the 
observations. In particular, the agree ment improves whe n 
taking into account a clumping factor iGranato et al .12004) . 
which allows the gas to be efficiently converted into stars also 
in very massive dark haloes with SFR> 1000 Mq yr~^. 

In this way, narrow associated QSO absorption systems 
proved to be extremely useful in the study of the QSO en- 
vironment, in particular when there is evidence of their in- 
trinsicness. They can be used as estimators of the chemical 
abundances in high redshift spheroidal galaxies which are 
not easily determined otherwise. The probed gas will prob- 
ably be ejected from the galaxy due to the QSO feedback, 
thus we are also observing a potential source of enrichment 
of the intergalactic medium at high redshift. 

In order to obtain a deeper insight in the evolution of 
QSO host-galaxies and environments it is essential to enlarge 
the data sample. In particular, obtaining high signal-to-noise 
spectra in the UV to reliably measure the doubly-ionised 
lines of C and N and increasing the redshift range especially 
at la rge values. Indeed, th e five 2 ~ 4 AALs analised up to 
now JSavaglio et alJll997ft seems to indicate a slightly lower 
average metallicity, [C/H] ~ —0.5, than for the bulk of the 
sample at redshift z ^ 2 — 2.5. More data will be fundamen- 
tal to verify the observed correlations and to constrain the 
predictions of theoretical models. 
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Table 9. Results of the analysis of the AALs in our sample. Av is the velocity separation between the absorption system and the QSO 
emission redshift, Ar(Hi) is the neutral hidrogen column density of the studied component and U is the ionisation parameter. We adopted 
the standard definition: [X/Y] = log (X/Y) - log (X/Y)q and the solar abundances determined by Holweger (200f). Uncertainties on 
the relative abundance determinations are f a. 



QSO UM680 UM68f HEff58-f843 Q2343+f 232(1) Q2343+1232(2) Q2343+1232(3) Q0453-423 PKS0329-255 



^abs 


2.123 


2.122 


2.4426 


2.5688 


2.5695 


2.5699 


2.6362 


2.7089 


Av (km s^^) 


1990 


-10 


900 


-1668 


-1727 


-1760 


2040 


-1940 


logW(Hi) (cm-2) 


15.2 ± 0.6-" 


15.7 + 0.2 


14.22 + 0.02 


13.42 + 0.01 


13.65 + 0.02 


13.99 + 0.01 


14.76 + 0.01 


14.01 + 0.05 


log[/ 


-2.26 ± 0.03 


-2.42 + 0.04 


-0.47 + 0.03 


-1.56 + 0.06 


-1.76 + 0.04 


-1.69 + 0.04 


-1.63 + 0.03 


-0.63 + 0.05 


[C/H] 


+ 1.4 ± 0.6 


+0.0 + 0.2 


+0.30 ± 0.06 


+0.46 + 0.08 


+0.37 + 0.05 


+0.50 + 0.09 


-0.80 + 0.02 


-0.04 + 0.05 


[N/H] 


+ 1.4 ± 0.6 


+0.0 + 0.2 


+0.38 ± 0.05 


+0.9 + 0.1 


+0.84 + 0.08 


+0.9 + 0.1 


-1.40 + 0.07 


-0.46 + 0.05 


[N/C] 


+0.0 ± 0.2 


+0.0 + 0.1 


+0.08 + 0.06 


+0.5 + 0.1 


+0.47 + 0.09 


+0.4 + 0.1 


-0.60 + 0.07 


-0.42 ± 0.02 


[O/C] 
















+0.2 + 0.1 


[Mg/C] 


+0.2 ± 0.2 


+0.6 + 0.1 














[Al/C] 


-1.0 ± 0.2 


+0.50 + 0.08 














[Si/C] 


-0.45 + 0.2 


+0.5 + 0.1 








+0.3 + 0.1 


+0.70 + 0.04 





" See the discussion about this system in Section 3.2 for an explanation on how this column density and its error were determined. 



